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Lattice QCD calculations of two-nucleon systems are used to isolate the short-distance two-body 
electromagnetic contributions to the radiative capture process np —>■ d-y, and the photo-disintegration 
processes —>■ np. In nuclear potential models, such contributions are described by phenomeno¬ 

logical meson-exchange currents, while in the present work, they are determined directly from the 
quark and gluon interactions of QCD. Calculations of neutron-proton energy levels in multiple 
background magnetic fields are performed at two values of the quark masses, corresponding to pion 
masses of mT, ~ 450 and 806 MeV, and are combined with pionless nuclear effective field theory to 
determine these low-energy inelastic processes. Extrapolating to the physical pion mass, a cross sec¬ 
tion of > d'y) = 332.4( 'f’fj ) mb is obtained at an incident neutron speed of u = 2, 200 m/s, 

consistent with the experimental value of a^^^^{np —>■ d'y) = 334.2(0.5) mb. 

PACS numbers: 11.15.Ha, 12.38.Gc, 13.40.Gp 


The radiative capture process, np —)■ d'y, plays a crit¬ 
ical role in big-bang nucleosynthesis (BBN) as it is the 
starting point for the chain of reactions that form most 
of the light nuclei in the cosmos. Studies of radia¬ 
tive capture m, and the inverse processes of deuteron 
electro- and photo-disintegration, 'y^*'>d —>■ np mzi, have 
constrained these cross-sections and have also provided 
critical insights into the interactions between nucleons 
and photons. They conclusively show the importance of 
non-nucleonic degrees of freedom in nuclei, which arise 
from meson-exchange currents (MECs) in the context 
of nuclear potential models mm- Nevertheless, in the 
energy range relevant for BBN, experimental investiga¬ 
tions are challenging [10]. For the analogous weak in¬ 
teractions of multi-nucleon systems, considerably less is 
known from experiment but these processes are equally 
important. The weak two-nucleon interactions currently 
contribute the largest uncertainty in calculations of the 
rate for proton-proton fusion in the Sun EHn], and in 
neutrino-disintegration of the deuteron [18], which is a 
critical process needed to disentangle solar neutrino os¬ 
cillations. Given the phenomenological importance of 
electroweak interactions in light nuclei, direct determi¬ 
nations from the underlying theory of strong interaction, 
quantum chromodynamics (QCD), are fundamental to 
future theoretical progress. Such determinations are also 
of significant phenomenological importance for calibrat¬ 
ing long-baseline neutrino experiments and for investiga¬ 


tions of double beta decay in nuclei. In this Letter, we 
take the initial steps towards meeting this challenge and 
present the first lattice QCD (LQCD) calculations of the 
np —>■ d'y process. The results are in good agreement with 
experiment and show that QCD calculations of the less 
well-determined electroweak processes involving light nu¬ 
clei are within reach. Similarly, the present calculations 
open the way for QCD studies of light nuclear matrix ele¬ 
ments of scalar [19] (and other) currents relevant for dark 
matter direct detection experiments and other searches 
for physics beyond the Standard Model. 

The low-energy cross section for np —> d'y is conve¬ 
niently written as a multipole expansion in the electro¬ 
magnetic (EM) field piTl 1^ . 

+ ( 1 ) 

where Xmi is the Ml amplitude, 70 is the binding mo¬ 
mentum of the deuteron, M is the mass of the nucleon, 
and p is the momentum of each incoming nucleon in the 
center-of-mass frame. The ellipsis denotes the contribu¬ 
tion from El and higher-order multipoles (higher multi¬ 
poles can be included systematically and improve the re¬ 
liability of the description [52] , but are not relevant at the 
level of precision of the present work). In a pionless effec¬ 
tive field theory expansion [231 - 125] . employing dibaryon 
fields to resum effective range contributions [HHZ], the 
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leading-order (LO) and next-to-leading order (NLO) con¬ 
tributions lead to the M1 amplitude [27l [28] 


Xmi = 


-ir + -*lpl 
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where ki = (np — Kn) /2 is the isovector nucleon mag¬ 
netic moment, Zd = l/-\/l ~ 7 o ’'3 is the square-root of 
the residue of the deuteron propagator at the pole with 
ra the effective range in the channel, and ai, ri are the 
scattering length and effective range in the ^S'o channel. 
The quantity h = h — y/rir^Ki encapsulates the short- 
distance two-nucleon interactions through /i, but also de¬ 
pends on Ki- It is well established that gauge-invariant 
EM two-nucleon interactions (and direct photon-pion 
couplings in pionful effective held theories) [HI [THl [Hj- 
[M1 must be included in order to determine radia¬ 

tive capture and breakup cross-sections to a precision of 
better than ^ 10%. 

The only quantity in Eqs. 0 and 0 that is not deter¬ 
mined by kinematics, single-nucleon properties or scat¬ 
tering parameters, is h. In this work, we use LQCD 
to calculate this quantity by determining the energies of 
neutron-proton systems in background magnetic helds. 
A magnetic held mixes the Iz = jz = 0 n-p states in the 
and ^Si-^Di channels, providing sensitivity to the 
EM interactions. The deuteron and dineutron ground 
states are nearly degenerate at both pion masses used in 
the present calculation [33|, and the two-nucleon sector 
exhibits an approximate spin-havor SU(4) symmetry (as 
predicted by the large-A), limit of QCD [31]). In this case, 
it can be shown |35j that the energy diherence between 
the two eigenstates depends upon li as 


AEss,,^Soi^) = ‘^{^i+loZjh) ^|B|+0(|Bn, (3) 


where B is the background magnetic held. It is con¬ 
venient to focus on the combination Li — that 

characterizes the two-nucleon contributions. 

Our LQCD calculations were performed on two en¬ 
sembles of gauge-held conhgurations generated with a 
clover-improved fermion action |36j and a Liischer-Weisz 
gauge action m- The hrst ensemble had Nf = 3 de¬ 
generate light-quark havors with masses tuned to the 
physical strange quark mass, producing a pion of mass 
m,r 806 MeV and used a volume oi xT = 32^ x 48. 
The second ensemble had Nf = 2 + 1 havors with the 
same strange quark mass and degenerate up and down 
quarks with masses corresponding to a pion mass of 
450 MeV and a volume oi x T = 32^ x 96. 
Both ensembles had a gauge coupling of (3 = 6.1, cor¬ 
responding to a lattice spacing of a ^ 0.12 fm. Back¬ 
ground EM (t/Q(l)) gauge helds giving rise to uniform 
magnetic helds along the xa-axis were multiplied onto 


each QCD gauge held in each ensemble (separately for 
each quark havor), and these combined gauge helds were 
used to calculate up-, down-, and strange-quark propa¬ 
gators, which were then contracted to form the requi¬ 
site nuclear correlation functions using the techniques 
of Ref. [3H]- Calculations were performed on ^ 1,000 
gauge-held conhgurations at the SU(3) point and ~ 650 
conhgurations at the lighter pion mass, each taken at in¬ 
tervals of 10 hybrid Monte-Carlo trajectories. On each 
conhguration, quark propagators were generated from 48 
uniformly distributed Gaussian-smeared sources for each 
magnetic held. For further details of the production at 
the SU(3)-symmetric point, see Refs. [33l|39l|40] and in 
particular, Ref. [35]. Analogous methods were employed 
for the calculations using the lighter pion mass ensemble. 

Background EM helds have been used extensively to 
calculate electromagnetic properties of hadrons, such as 
the magnetic moments of the lowest-lying baryons gU- 
[49] and light nuclei [40] , and the polarizabilities of mesons 
and baryons gHl 1101 • The quark helds have electric 
charges Qu = +2/3 and Qd,s = —1/3 for the up-, down- 
and strange-quarks, respectively, and background mag¬ 
netic helds are required to be quantized m in order that 
the magnetic hux is uniform throughout the lattice. The 
link helds, Up^\x), associated with the background held 
are of the form 


= e 
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for quark havor g, where h is an integer. The uni¬ 
form magnetic held resulting from these links is e B = 
dirn/L^z, where e is the magnitude of the electric charge 
and z is a unit vector in the xa-direction. In physical 
units, the background magnetic helds used with these en¬ 
sembles of gauge conhgurations are e|B| ^ 0.05|h| GeV^. 
To optimize the re-use of light-quark propagators in the 
calculations, Uq{1) helds with h = 0,1, —2,4 were used. 
At the SU(3) symmetric point, additional calculations 
were performed with h = 3, —6,12. 

With three degenerate havors of light quarks, and a 
traceless electric-charge matrix, there are no contribu¬ 
tions from the magnetic held coupling to sea quarks at 
the SU(3) point at leading order in the electric charge. 
This is not the case for the ~ 450 MeV calculations 
because of havor SU(3) breaking. However, Li is an 
isovector quantity in which sea quark contributions can¬ 
cel (the up and down sea quarks used in this work are 
degenerate) so it is correctly determined by the present 
calculations. 

In this work, we focus on the Iz = Jz = 0 coupled- 
channel neutron-proton systems. Our analysis follows 
that of Ref. [35] which presents results on the "^806 
MeV ensemble, and we direct the reader to that work 
for more detail regarding the interpolating operators and 
statistical analysis methods that are used. A matrix of 
correlation functions generated from source and sink op- 
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FIG. 1: The double ratios of the two principal correlators 
are shown for ~ 450 MeV for the three magnetic field 
strengths. The bands correspond to the single-exponential fits 
to the correlator and the associated statistical uncertainty. 



FIG. 2: LQCD calculations of the energy-splittings between 
the two lowest-lying eigenstates, with the single-nucleon con¬ 
tributions removed, as a function of h, along with the asso¬ 
ciated fits. The lower (blue) set of points correspond to the 
niTr ~ 450 MeV ensemble and the upper (green) points to 
m,r ~ 806 MeV. The slope of the sets of points is propor¬ 
tional to Li at the appropriate pion mass. 


erators associated with and ^ 5*0 Iz = jz = 0 interpo¬ 
lating operators 


C(t;B) 




is diagonalized to yield “principal correlators”, A±(t;B), 
corresponding to the eigenstates of the coupled sys¬ 
tem. In all cases, the principal correlators exhibit single¬ 
exponential behavior at times where statistical uncertain¬ 
ties are manageable. To highlight the difference arising 
from purely two-body effects, a ratio of ratios of the prin¬ 
cipal correlators to the appropriate single particle corre¬ 
lation functions is formed 




A+(t;B)C„,^(t;B)Cp,;(t;B) 
A_(t;B)C„,^(t;B)Cp.t(t;B) ’ 


FIG. 3: The results of LQGD calculations of Li (blue 

points). The blue (green) shaded regions point show the lin¬ 
ear (quadratic) in extrapolation of Li to the physical pion 
mass (dashed line) in natural nuclear magnetons (nNM). The 
vertical (red) line indicates the physical pion mass. 


where are the correlation functions corre¬ 

sponding to the different polarizations of the proton and 
neutron. For large time separations, 

SR^Su^Soit oo;B) A ( 7 ) 

where A is an overlap factor and the energy shift is 

SEas^iSo = ‘^E^Si,^So - [Ep,t ~ Rp,l] + [Rn.t - En,l] 

^ 2Li|eB|/M + 0(B2), (8) 

omitting the B dependence for clarity. Fig. shows 
the above ratios for the ttItt ~ 450 MeV ensemble for 
each magnetic field strength, along with correlated single¬ 
exponential fits to the time dependence and their statisti¬ 
cal uncertainties. The energies extracted from these fits 
depend on |B|, with ‘2-^Li being the coefficient of the 
linear term. Fig. shows the extracted energy shifts for 
both the m-TT ~ 450 MeV and 806 MeV ensembles. The 
figure also shows the envelopes of a large range of poly¬ 
nomial fits to their magnetic field dependence. Ref. [3S] 
presents the WTr ~ 806 MeV correlation functions in de¬ 
tail, and has a complete discussion of the fitting methods 
used in the analysis for both sets of pion masses. 

The extracted values of Li are shown in Fig.j^for both 
sets of quark masses. The functional dependence of Li 
on the light-quark masses is not known. However, the 
deuteron and dineutron remain relatively near threshold 
over a large range of quark masses [33l [52H^ , and the 
magnetic moments of the nucleons are essentially inde¬ 
pendent of the quark masses when expressed in units of 
natural nuclear magnetons m, so it is plausible that 
Li also varies only slowly with the pion mass. Indeed, 
there is only a small difference in the value of Li at 
~ 806 MeV and at ~ 450 MeV. In order to 
connect to the physical point, we extrapolate both lin¬ 
early and quadratically in the pion mass by resampling 
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the probability distribution functions of Li determined 
by the field-strength dependence fits at each pion mass. 
The two forms of extrapolation yield consistent values 
at the physical point, with the central value and uncer¬ 
tainties determined from the 0.17, 0.50 and 0.83 quan¬ 
tiles of the combination of the two projected probabil¬ 
ity distribution functions. After this extrapolation, the 
value = 0.285( ^eo ) nNM is found at the physical 
pion mass, where the uncertainty incorporates statisti¬ 
cal uncertainties, correlator fitting uncertainties, field- 
strength dependence fitting uncertainties, and the uncer¬ 
tainties in the mass extrapolation. This leads to a value 
^iqcd _ +16 ^ Future calculations with lighter 

quark masses will reduce both the statistical and system¬ 
atic uncertainties associated with Li. 

The cross section for np —>■ dy has been precisely mea¬ 
sured in experiments at an incident neutron speed of 
u = 2, 200 m/s [1]. Using the expressions in Eqs. Q and 
([^, the experimentally determined deuteron binding en¬ 
ergy and ^S'o scattering parameters, the experimentally 
determined nucleon isovector magnetic moment, and the 
above extrapolated LQCD value of leads to a cross 
section at u = 2,200 m/s of 

= 332.4( ) mb , (9) 

which is consistent with the experimental value of = 
334.2(0.5) mb [1] within uncertainties (see also. Ref. [56]). 
As in the phenomenological determination, the two-body 
contributions are 0(10%). At the quark masses where 
the lattice calculations are performed, the cross-sections 
are considerably smaller than at the physical point, pri¬ 
marily because the deuteron binding energy is larger. At 
806 MeV, the scattering parameters, binding en¬ 
ergy and magnetic moments have been determined previ¬ 
ously [MiiMiiin] and we can predict the scattering cross 
section using only lattice QCD inputs, with a median 
value ~ 5 mb at u = 2, 200 m/s.^ 

Summary: Lattice QCD calculations have been used 
to determine the short-distance two-nucleon interactions 
with the electromagnetic field (meson-exchange currents 
in the context of nuclear potential models) that make sig¬ 
nificant contributions to the low-energy cross-sections for 
np —>■ d'j and > np. This was facilitated by the pio- 

nless effective field theory which provides a clean separa¬ 
tion of long-distance and short-distance effects along with 
a concise analytic expression for the near-threshold cross 
sections. A (naive) extrapolation of the LQCD results 
to the physical pion mass is in agreement with the ex¬ 
perimental determinations of the np —>■ dy cross-section, 
within the uncertainties of the calculation and of the ex¬ 
periment. Calculations were performed at a single lattice 


^ Propagation of the uncertainties in the required inputs leads to 
a highly non-Gaussian distribution of 1351 . 


spacing and volume, introducing systematic uncertainties 
in Li that are expected to be small in comparison to our 
other uncertainties, ©(o^Aqqq, 6“'’'“'^) < 4%. A 

more complete study, and a reduction of the uncertainties 
of this cross-section will require additional calculations at 
smaller lattice spacings and larger volumes, along with 
calculations at smaller quark masses. 

The present calculation demonstrates the power of lat¬ 
tice QCD methods to address complex processes of im¬ 
portance to nuclear physics directly from the Standard 
Model. The methods that are used are equally applica¬ 
ble to weak processes such as pp —)■ de^v^ vd _ppe+, 
vd —>■ ^d, and vd —>■ vnp, as well as to higher-body tran¬ 
sitions. Background field techniques will also enable the 
extraction of nuclear matrix elements of other currents 
relevant for searches for physics beyond the Standard 
Model. Extensions of our studies to larger systems are 
currently under consideration, and calculations in back¬ 
ground axial-vector fields necessary to address weak in¬ 
teraction processes are under way. As this technique has 
successfully recovered the short-distance contributions to 
np —>■ dy, it also seems likely that it can be generalized 
to the calculation of parity-violating observables in this 
process resulting from weak interactions, or from physics 
beyond the Standard Model (see Ref. m for a review). 
Einally, the present work reinforces the utility of com¬ 
bining lattice QCD calculations with low-energy effective 
field theories describing multi-nucleon systems [58] . 
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